Complete relativistic second-order dissipative hydrodynamics from the entropy 

principle 
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We present a new derivation of relativistic dissipative hydrodynamic equations, which invokes the 
second law of thermodynamics for the entropy four-current expressed in terms of the single-particle 
phase-space distribution function obtained from Grad's 14-moment approximation. This derivation 
is complete in the sense that all the second-order transport coefficients are uniquely determined 
within a single theoretical framework. In particular, this removes the long-standing ambiguity in 
the relaxation time for bulk viscosity thereby eliminating one of the uncertainties in the extraction 
of the shear viscosity to entropy density ratio from confrontation with the anisotropic flow data 
in relativistic heavy-ion collisions. We find that in the one-dimensional scaling expansion, these 
transport coefficients prevent the occurrence of cavitation even for rather large values of the bulk 
viscosity estimated in lattice QCD. 
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Relativistic fluid dynamics has been quite successful in 
explaining the various collective phenomena observed in 
cosmology, astrophysics and the physics of high-energy 
heavy-ion collisions. The earliest theories of relativis- 
tic dissipative hydrodynamics by Eckart and Landau- 
Lifshitz ^2j] were based on the assumption that the en- 
tropy four-current is first order in dissipative quantities, 
which led to parabolic differential equations that suf- 
fered from acausality. The second-order Israel-Stewart 
(IS) theory Q with the entropy current quadratic in dis- 
sipative quantities led to hyperbolic equations and thus 
restored causality. 

Application of the second-order viscous hydrody- 
namics to high-energy heavy-ion collisions has evoked 
widespread interest ever since a surprisingly small value 
for the shear viscosity to entropy density ratio rj/s was 
estimated from the analysis of the elliptic flow data [3| ■ 
Indeed the estimated r]/s was close to the conjectured 
lower bound ry/sjKSS = l/47r This led to the claim 

that the quark-gluon plasma (QGP) formed at the Rela- 
tivistic Heavy-Ion Collider (RHIC) was the most perfect 
fluid ever observed. A precise estimate of rj/s is vital to 
the understanding of the properties of the QCD matter. 

In this Communication, we provide a solution to one 
of the major uncertainties that hinders an accurate ex- 
traction of the viscous corrections to the ideal fluid be- 
havior, namely the inadequate knowledge of the second- 
order transport coefficients. In the standard derivation 
of second-order evolution equations for dissipative quan- 
tities from the requirement of positive divergence of the 
entropy four-current, the most general algebraic form of 
the entropy current is parameterized in terms of unknown 
thermodynamic coefficients Q. These coefhcients which 
are related to relaxation times and coupling lengths of 
the shear and bulk pressures and heat current, however, 
remain undetermined within the framework of thermo- 
dynamics alone While kinetic theory for massless 



particles and strongly coupled A/" — 4 supersymmetric 
Yang-Mills theory Q predict different shear relaxation 
times = 3/2ttT and (2 — ln2)/27rT, respectively, for 
rj/s = l/47r, the bulk relaxation time rn remains com- 
pletely ambiguous. Hence ad hoc choices have been made 
for the value of rn in hydrodynamic studies lol- 13|. 

Lattice QCD studies for gluonic plasma in fact pre- 
dict large values of bulk viscosity to entropy density ra- 
tio, C/s, of about (6-25) ?7/s|kss near the QCD phase- 
transition temperature Tc 1J| . This would translate into 



large values of the bulk pressure and bulk relaxation time, 
and may affect the evolution of the system significantly 
12[. Further, the large bulk pressure could result 
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in a negative longitudinal pressure leading to mechani- 
cal instabilities (cavitation) whereby the fluid breaks up 
into droplets [13, [H, [3 ■ Thus the theoretical uncertain- 
ties arising from the absence of reliable estimates for the 
second-order transport coefficients should be eliminated 
for a proper understanding of the system evolution. 

We present here a formal derivation of the dissipa- 
tive hydrodynamic equations where all the second-order 
transport coefficients get determined uniquely within a 
single theoretical framework. This is achieved by in- 
voking the second law of thermodynamics for the gen- 
eralized entropy four-current expressed in terms of the 
phase-space distribution function given by Grad's 14- 
moment approximation. Significance of these coefficients 
is demonstrated in one-dimensional scaling expansion of 
the viscous medium. 

Hydrodynamic evolution of a medium is governed by 
the conservation equations for the energy-momentum 
tensor and particle current 17 [ 



T"'' ^ J dp p^'p^U + f)= eu'^u'' - (P + H)A'"' + ir"" , 
= I dp p^(/ - /) = nu^ + n^, (1) 
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where dp = (7dp/[(27r)^A/p^ 



], g and m being the de- Boltzmann's H- function: 



generacy factor and particle rest mass, is the particle 
four-momentum, / = f{x,p) is the phase-space distribu- 
tion function for particles and / for antiparticles. The 
above integral expressions assume the system to be di- 
lute so that the effects of interaction are small 17[. In 



the above tensor decompositions, e, P, n are respectively 
energy density, pressure, net number density, and the 
dissipative quantities are the bulk viscous pressure (11), 
shear stress tensor (tt^'^) and particle diffusion current 
(n'^). Here A^'^ = g'^'^ — u'^u'^ is the projection operator 
on the three-space orthogonal to the hydrodynamic four- 
velocity u'^ defined in the Landau frame: T^'^Ui, — eu^^. 

Energy-momentum conservation, d^T^^ = and cur- 
rent conservation, d/^iN'^ = yield the fundamental evo- 
lution equations for e, u'^ and n. 



De + ie + P + n)a^7i^ - ^^''V(^w,) = 0, 
(e + P + n)Du" - V"(P n) + A^d^TT"" = 0, 

Dn + ndf.u'' + d^n'' = 0. (2) 

We use the standard notation ^("S^) = (A^B'^ + 
A^B")/2, D = and V" = A'^"9^. Even if the 

equation of state is given, the system of Eqs. ^ is not 
closed unless the evolution equations for the dissipative 
quantities H, tt^'', are specified. 

Traditionally the dissipative equations have been ob- 
tained by invoking the second law of thermodynamics, 
viz., d^S^ > 0, where the entropy four-current is 
given by |,[3,i 



(3) 



Here l3 — 1 /T is the inverse temperature, /i is the chemi- 
cal potential, a — /3/i, and is a function of deviations 
from local equilibrium. The second equality is obtained 
by using the definition of the equilibrium entropy density 
s = P[e + P ~ fin) and Taylor-expanding Q'' to second 
order in dissipative fluxes. In this expansion, /3i(e, n) > 
and ai(e, n) > are the thermodynamic coefficients cor- 
responding to pure and mixed terms. These coefficients 
can be obtained within the kinetic theory approach such 
as the IS theory 3]. However, it is important to note 
that they cannot be determined solely from thermody- 
namics using Eq. (|3]) and as a consequence the evolution 
equations remain incomplete. 

In contrast to the above approach, our starting point 
for the derivation of the dissipative evolution equations 
is the entropy four-current expression generalized from 



ql-l 



I dpp^ [/(ln/-l) + (/-^/)] 



5,^=±i = - I dpp^ ( / In / + rf In / ) + (/ ^ /) ,(4) 



where / = 1 — r/ and r = 1, —1,0 for Fermi, Bose, and 
Boltzmann gas, respectively. The divergence of ^.^ 
leads to 



9^5^ = - 



dpp^ {d,f)Hf/f) + {f f) . (5) 



For small departures from equilibrium, / and / can be 
written as f = fo + 6f and f ^ fa + Sf. The equilibrium 
distribution functions are defined as /o — [exp(/3u • p — 
a) + r]~^ and Jo = [exp(/3u-p-|-a)-|-r]"^, where /3 = 1/T 
and a = fi/T are obtained from the equilibrium matching 
conditions n = uq and e = eg- 

To proceed further, we take recourse to Grad's 14- 
moment approximation [l8l | for the single particle dis- 
tribution in orthogonal basis 13, 2^ 



f = h+fo.M, = Ann+A„nQp"-HA^7rQ/3p>^, (6) 

and similarly for /. The coefficients (An,A„,A7r) are as- 
sumed to be independent of four-momentum p^ and are 
functions of (e,a,/3). From Eqs. ^ and ([6]), we get 



dp 



id J) In 



In 1 



1 - rfo<P 



+ (/ ^ /, fo ^ /o) 



(7) 



The </)-independent terms on the right vanish due to 
energy-momentum and current conservation equations. 
To obtain second-order evolution equations for dissipa- 
tive quantities, one should consider 5"^ up to the same 
order. Hence d^S^ necessarily becomes third-order. Ex- 
panding the 0-dependent terms in Eq. ([7]) and retaining 
all terms up to third order in gradients (where 4> is linear 
in dissipative quantities), we get 



dp [{0 (a^/o) - 02(/o - i/2)(a^/o) 

+ 0'a^(/o/o) + cpfofoid^cj,)} + (fo ^ fo)] ■ (8) 

The various integrals in the above equation can be de- 
composed into hydrodynamic tensor degrees of freedom 
via the definitions: 



dpp/'i...p^"(/o±/o)=/±u''i...^M. 



perms) -|- 



(9) 



where 'perms' denotes all non-trivial permutations of 
the Lorentz indices. We similarly define j^i''^ "^" ^nd 
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j;j^MiM2 - A»n -^j-^gj-g ^}jg rnomentum integrals are weighted 

with /o/o ± (/o fo) and /o/o ± (/o /o), and are 
tensor decomposed with coefheients J^^ and K^^, respec- 
tively. All these coefficients can be obtained by suitable 
contractions of the integrals and are related to each other 
by 



2K^ 



1 



and also satisfy the differential relations 

^-^nq '^nq d/S nq ' dcX 



(10) 



d/3 



n — 1 , g 



da 



(11) 



With the help of these relations and Grad's 14-monient 
approximation, Eq. ([5]) reduces to 



-/?n 



+l3on +l3unUe + aoV.n^ + ^a„nrv^^ 



(12) 



where a^, axY, /3_xx are known functions of /?, a 
and the integral coefficients JjJ"^ and i^^^. Two new 

parameters "0 and x with — 1 — and X — 1 ~ X are 
introduced to 'share' the contributions stemming from 
the cross terms of 11 and Tr**" with n^. 

The second law of thermodynamics, d^S'^ > 0, is guar- 
anteed to be satisfied if we impose linear relationships be- 
tween thermodynamical fluxes and extended thermody- 
namic forces, leading to the following evolution equations 
for bulk, charge current and shear 



n 



-I- ipUniinfj.u^ + -0annnp V'a 



(13) 



+ XaTrnTT^^V'a + Xan7r<'u'' 



2V 



(14) 



(15) 



^(m> = /^i^A" and = A^^B"'' represent space-like 

and traceless symmetric projections respectively, both 
orthogonal to m'', where A'^''^ = [A^A^ + A^A^ - 
(2/3)A^''Ac,/3]/2. It may be noted that although the 
forms of the Eqs. P^ -([T5 |) are the same as in the stan- 
dard Israel- Stewart theory [1,01, all the transport coeffi- 
cients are explicitly determined in the present derivation: 



ao XnXnJ}i/f3, ol\ 



J3+//3, /32 = 2Ay5+//3, 
= -2A„A„ J4+2//3. (16) 



As a consequence, the relaxation times defined as, 

rn=C/3o, T„ = A/3i, t^=2^^2. (17) 



can be obtained directly. With An 



I/J21, A^ = l/(2J4+2), n = /rn, e = /o+n, and P = -1+ 



-'42^ — -"10' ~ -'20, 

the expressions for /3i,cko,CKi simplify to 



21, 



l3i = {e + P)/n'^, aa = ai = l/n. (18) 

For a classical Boltzmann gas {fo ~ 1), the coefficients 
/3o and /32 take the simple forms 

/3o = l/P, /32 = 3/(e + P) + m^l3^P/[2{e + Pf]. (19) 



Equations (|T3 l) -([T5 |) in conjunction with the second-order 
transport coefficients ([T5)) and constitute one of the 
main results in the present work. These coefficients are 
obtained consistently within the same theoretical frame- 
work. In contrast, in the standard derivation from en- 
tropy principles [Z] , the transport coefficients have to be 
estimated from an alternate theory. For instance, in the 
IS derivation based on kinetic theory, these involve com- 
plicated expressions which in the photon limit (m/3 — >■ 0) 



reduce to 21 



= 216/(to'*/3'*P), /3^^ = 3/4P. (20) 

An alternate derivation from kinetic theory (KT) using 
directly the definition of dissipative currents yields [llj 



KT 



1 



(e + P)--(6-3P) 



KT 



AP 
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{e-W) 
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where Cg is the speed of sound and (• • • ) = J dp{- 
A field-theoretical (FT) approach gives [23| 



(21) 
•)/o. 



/3o 



FT 



{e + P)--{e-'SP) 
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/3r=l/[2(3-a)P], 



(22) 



with the coefficients of charge conductivity, bulk and 
shear viscosity, viz. A, C,?? > 0. The notations. 



where a = 2 for charged scalar bosons and a = 3 for 
fermions. We find that our expression for /?2 (Eq. ((T^ ) 



T (GeV) 

FIG. 1; (Color online) Temperature dependence of bulk and 
shear relaxation times. Inset shows (/s (see text) and rj/s = 
l/47r. 



in the masslcss limit, agrees with the IS result (Eq. (I^H)) ) 
and also with those obtained in Refs. [1, [l^l. Thus the 
shear relaxation times (Eq. (|17p ) obtained here and in 
these studies are also identical. As Po in Eqs. ([201)- (EH) 
diverge in the massless limit, so does the bulk relaxation 
time Tn (Eq. pT]) ). thereby stopping the evolution of 
the bulk pressure. It is important to note that /3o in 
Eq. ()19|) and hence Tn in the present calculation remain 
finite in this limit. A detailed comparison of IS, KT and 
FT results can be found in 
■0 and X occurring in Eq. E 
in 



2J|. The two parameters 



remain undetermined as 
however, these do not contribute to the scaling 
expansion. 

To demonstrate the numerical significance of the new 
coefficients derived here, we consider the evolution equa- 
tions in the boost-invariant Bjorken hydrodynamics at 
vanishing net baryon number density ,25]. In terms 
of the coordinates (r, a;,?/,ry) where r = \/t'^ — and 
r\ — tanh^^(z/i), the initial four- velocity becomes = 
(1,0,0,0). For this scenario = and the evolution 
equations for e, tt = — t^tt'''' and 11 reduce to 



dr 


1 

T 


(e + p + n 


dix 


4?/ 


4t^ 








r 


3r 



n. 



(23) 
(24) 
(25) 



Noting that = 1/ P ^ (i^ = ^/ {e + P) and s = {e + P)/T, 
the relaxation times defined in Eq. (|17p reduce to 



p fc 



PT 



6 (rj 
s / I \s 



(26) 



We have used the state-of-the-art equation of state [26j , 
which is based on a recent lattice QCD result [27|. For 
C/s at T > Tc ~ 184 MeV, we use the parametrized 
form [13] of the lattice QCD results of Meyer [l4?] which 
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FIG. 2: (Golor online) (a) Time evolution of shear stress in 
the absence of bulk (11 = 0) and magnitude of bulk stress 
for Tn = C/P s-nd '''n = tv. The arrow indicates the time 
when Tc is reached, (b) Temperature dependence of pressure 
anisotropy, Pl/Pt, for these three cases. The results are for 
initial T = 310 MeV, to = 0.5 fm/c and ri/s = l/47r. The 
evolution is stopped when Pl vanishes. 



suggest a peak near Tc. At T < the sharp drop in 
^/s reflects its extremely small value found in the hadron 
resonance gas model [28]; see inset of Fig. [1] For the rj/s 
ratio, we use the minimal KSS bound [6| value of l/47r. 

In the absence of any reliable prediction for the bulk re- 
laxation time Tn , it has been customary to keep it fixed or 
set it equal to the shear relaxation time Tjt or parametrize 
it in such a way that it captures critical slowing-down of 
the medium near Tc due to growing correlation lengths 
13" 3]. Since C./s has a peak near the phase transition, 
the Tn obtained here (Eq. (j26|) ) and shown in Fig. [1] nat- 
urally captures the phenomenon of critical slowing-down. 

The evolution equations ([23|) - ([25|l are solved simulta- 
neously with an initial temperature Tq — 310 MeV jisj 
and initial time tq = 0.5 fm/c typical for the RHIC en- 
ergy scan. We take initial values for bulk stress and shear 
stress, n = TT = GeV/fm^ which corresponds to an 
isotropic initial pressure configuration. 

Figure [21^ a) shows time evolution of the shear pressure 
TT and the magnitude of the bulk pressure 11. At early 
times T < 2 fm/c or equivalently at T > 1.2Tc, shear 
dominates bulk. This implies that eccentricity-driven el- 
liptic flow which develops early in th e sy stem would be 
controlled more by the shear pressure [ij] . At later times 
(when T ^ Tc), the large value of (/s makes the bulk 
pressure dominant. This leads to sizeable entropy gen- 
eration (Eq. (|12p ) and consequently enhanced particle 
production. 

Figure [2fa) also compares the 11 evolution for bulk re- 
laxation time, Tn, calculated from Eq. ([26]) (solid line) 
and Tn = t^ (dashed line). At early times, the larger 
value of Tn in the latter case (see Fig. [1]) results in a 
relatively smaller growth of jll] as evident from Eq. 
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Near Tc, the rapid increase in C/s causes |n| to increase. 
Subsequently the longitudinal pressure Pr = (P + n — tt) 
vanishes leading to cavitation [l^, [13, Eal- con- 
trast, with our Tjj, this rise in (/s is overcompensated by 
a faster increase in rn thereby slowing down the evolu- 
tion of n. This behavior prevents the onset of cavitation 
and guarantees the applicability of hydrodynamics with 
bulk and shear up to temperatures well below Tc into the 
hadronic phase. Furthermore, this slowing down of the 
medium followed by its rapid expansion, has the right 
trend to explain the identical-pion correlation measure- 
ments (Hanbury Brown- Twiss puzzle) [2^, . 

The absence of cavitation in our calculation is clearly 
evident in Fig. EJb) which shows the variation of pressure 
anisotropy, Pl/Pt = (P+n-7r)/(P+n+7r/2), with tem- 
perature. Near Tc, the longitudinal pressure Pl vanishes 
if one assumes Tn = tv (dashed line) leading to cavita- 
tion, whereas it is found to be positive for all tempera- 
tures with Tn derived here (solid line). In fact, we have 
found that in the latter case, cavitation is completely 
avoided for the entire range of (/s values (0.5 < (/s < 2.0 
near Tc) estimated in lattice QCD T?]. The sizeable dif- 
ference between the 11 = case (dot-dashed line) and the 
Tn = C/P case (solid line) clearly underscores the impor- 
tance of bulk pressure near Tc, which can have significant 
implications for the elliptic fiow V2 III thus affecting the 
extraction of 77/s. Further, the large bulk pressure when 
incorporated in the freezeout prescription could also af- 
fect the final particle abundances and spectra. 

We have also found that the evolution of 11 is insensi- 
tive to the choice of initial conditions such as n(ro) = 
and the Navier-Stokes value — C(ro)/ro. This is due to 
very small Tn at early times (or higher temperatures) 
which causes 11 to quickly lose the memory of its initial 
condition and to relax to the same value at t ^ 1 fm/c. 

To summarize, we have presented a new derivation of 
the relativistic dissipative hydrodynamic equations from 
entropy considerations. We arrive at the same form of 
dissipative evolution equations as in the standard deriva- 
tion but with all second-order transport coefficients such 
as the relaxation times and the entropy flux coefficients 
determined consistently within the same framework. We 
find that in the Bjorken scenario, although the bulk pres- 
sure can be large, the relaxation time derived here pre- 
vents the onset of cavitation due to the critical slowing 
down of bulk evolution near Tc. 
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